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We study the deterministic dynamics of non-interacting classical gas particles confined to a one-dimensional 
box as a pedagogical toy model for the relaxation of the Boltzmann distribution towards equilibrium. Hard 
container walls alone induce a uniform distribution of the gas particles at large times. For the relaxation 
of the velocity distribution we model the dynamical walls by independent scatterers. The Markov property 
guarantees a stationary but not necessarily thermal velocity distribution for the gas particles at large times. 
We identify the conditions for physical walls where the stationary velocity distribution is the Maxwell dis- 
tribution. For our numerical simulation we represent the wall particles by independent harmonic oscillators. 
The corresponding dynamical map for oscillators with a fixed phase (Fermi-Ulam accelerator) is chaotic 
for mesoscopic box dimensions. 
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1 Introduction 

Classical mechanics and electrodynamics, quantum mechanics and quantum field theory are based on 
a deterministic time evolution. For given initial conditions, the particle positions and velocities or the 
values of the electromagnetic and Schrodinger fields are known for all times. Therefore, a lecturer of an 
introductory course on statistical mechanics faces the demanding task to justify a probabilistic description 
of physical processes on macroscopic time and length scales. 

Qualitative arguments which emphasize the importance of the thermodynamic limit are helpful. Stu- 
dents can certainly understand intuitively that information gets lost by destructive interference of individual 
signals: a man who listens simultaneously to ten thousand discussing philosophers receives as much in- 
formation as if he was herding a flock of sheep. Nevertheless, it is always helpful to provide explicit 
calculations for simple examples where the student can explicitly see that the outcome of a measurement 
on a macroscopic deterministic system can be derived equivalently from the consideration of a probabilistic 
substitute system. 

In this work, we provide such a simple example. We consider a toy model of classical non-interacting 
point particles in a box. We assume that our initial conditions factorize in the three spatial coordinates. 
Therefore, without a loss of applicability to three dimensions, we focus on one spatial dimension. We 
investigate how an initially non-thermal Boltzmann distribution function relaxes to thermal equilibrium 
for large times. As we develop our model, we shall encounter the known requirements for a stochastic 
description of a deterministic system, e.g., the coupling of the gas particles to a reservoir. It is the purpose 
of this work to show step by step how the relaxation of independent point particles is accomplished, without 
referring to concepts of statistical mechanics such as the entropy or the Boltzmann StoBzahl Ansatz. For a 
thorough discussion of one-dimensional models which employ the Boltzmann equation, see Ref. (TJ. 

In Sect. [2] we specify the 'swimming-lane' model for the particle propagation and the prediction for the 
Boltzmann distribution from statistical mechanics. In Sect.|3]we consider hard walls where the distribution 
of particle velocities remains constant in time. We show analytically that the particles are homogeneously 
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distributed for long times. The relaxation of the particle velocities requires dynamical walls where the 
gas particles exchange energy with wall particles. In Sect. |4]we consider wall particles whose velocities 
are chosen from a probability distribution. They scatter elastically with the gas particles. We analytically 
determine necessary conditions for a 'physical wall' where the gas particles' velocity distribution relaxes 
to a Maxwell distribution. In Sect. [5] we investigate numerically the deterministic time evolution of the gas 
particles from their elastic scattering against the wall particles which we model as harmonic oscillators. A 
summary and conclusions, Sect. [6] closes our presentation. 

2 Model for classical particles in a box 

One of the simplest model systems for the relaxation to thermal equilibrium are ideal classical particles in 
a box. We consider a large number N ^> 1 of non-interacting point-particles gi of mass m which are fully 
characterized by their velocities vt and their positions xg (\xg\ < a) in a one-dimensional box with walls 
at | a; | = a. Since the particles do not interact with each other we can treat them separately. They move 
independently between the walls such as swimmers lap in their individual lanes (swimming-lane model). 

2. 1 Initial conditions for the box particles 

At time t = 0, we presume that all particles start at position x = in the middle of the box. Their starting 
velocities are distributed symmetrically, i.e., for each particle gi with velocity vt there is a particle g e with 
velocity vj = —V(. The probability distribution for the event to find a particle with velocity v is denoted 
by the positive, normalized function Po(v) = Po(—v) > 0, dvPo(v) = 1. The typical velocity of 
the particles at t = 0, v typ , obeys v^ yp = (w 2 )o = dvv 2 Po(v). We assume that all moments of the 
velocity distribution exist, (v n )o < oo. 

2.2 Predictions from statistical mechanics for large times 

In our very dilute gas, the interaction between the particles and the walls, not among the particles them- 
selves, is responsible for the thermal equilibrium at large times, as it is expressed, e.g., in the ideal gas 



Here, p is the pressure which the N particles exert onto the walls of the container of volume V = 2a; fee is 
the Boltzmann constant which permits us to express the energy content of the wall in terms of the familiar 
temperature T J2] . 

For our system of classical particles, classical statistical mechanics makes elaborate predictions be- 
yond ([TJ about the particles' properties in thermodynamic equilibrium. Let AN (x, v, t) be the number of 
particles which can be found in the mesoscopic position interval [x, x + Ax[ (Ax = a/y/W -C a) with 
velocities from the mesoscopic interval [v, v + Av[ (Av = v typ /\/N <C Vtyp) at time t. We write 



In the thermodynamical limit, TV — s- oo so that (Ax)(Av) — y 0, n(x, v, t) approaches the familiar Boltz- 
mann distribution function f(x, v, t) |3J. From the Boltzmann distribution function we obtain the distribu- 
tion functions for the particle positions W(x, t) and their velocities P(v, t) as 



law, 



pV = Nk B T . 



(1) 



AN(x, v, t) = Nn(x, v, t)(Ax)(Av) . 



(2) 




(3) 



At time t = 0, we start from 

f(x,v,t = 0) - 

so that W(x, t = 0) = S(x) and P(v, t = 0) = P (v). 



6{x)P (v) 



(4) 
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For large times, statistical mechanics predicts that the Boltzmann distribution approaches its thermal 
equilibrium f oq (x, v), 

f \x\^\ / 777 f inv^ \ 

f{x,v,t^oo) = Ux,v) = -e\l- T )Mv) ; f u (v) = ) J^ j ^c xp \ c —j , (5) 

where Q(x) is the Heaviside step-function. The particles are equally distributed in the box, W cq (x) = 
W(x, t — > oo) = (l/2a)0(l — |x|/a), and their velocity distribution is given by Maxwell's expression, 

P(v,t oo) = P cq (v) = f M (v). 



3 Relaxation of particle positions 

Thus far, we have not yet specified the properties of the walls. We start with hard walls so that the particles 
are perfectly reflected: Every particle gi with velocity v i which is reflected at the right wall (v £ > n- v' e = 
—vi) has its mirror particle ~g t with velocity vj = —vt which is reflected at the left wall (vj M- Vj = —vj = 
vt). Consequently, the velocity distribution of the particles does not change in time, P(v, t) = Pq(v). 

3.1 Free Boltzmann equation in the absence of walls 
If the walls were absent, all particles start at x = and propagate freely according to Newton's first law, 

xe(t) = v e t . (6) 

Consequently, the particles at position x with velocity v at time t arrive at position x + vAt with velocity 
v at time t + At. This free propagation in phase space implies 

n(x, v,t) =n(x + vAt,v,t + At) . (7) 

In the thermodynamic limit and for small At, the Taylor expansion of this equation implies the free Boltz- 
mann equation, 

M|M+3l^=0. (8) 
ot ox 

Its general solution is f(x, v, t) = g(x — vt, v). Together with the initial condition (0]) we find 

f fTee {x,v,t) = 6(x-vt)P (v) . (9) 

The velocity distribution does not change in time, P(v, t) = Po(v). The position distribution broadens in 
time, W(x,t) = P(x/t)/t. The particles which were initially concentrated at the origin, disperse into a 
cloud whose density goes to zero as a function of time. 



3.2 Free Boltzmann equation in the presence of hard walls 

In the presence of hard walls, we consider infinite replicas of the box on the whole x-axis. Box number n 
is equivalent to the interval [na — a,na + a[ (extended box scheme). In the extended box scheme, the 
particle propagate freely, see eq. ©. In order to determine their true position, we have to fold the extended 
box scheme back into the interval for tt, = (reduced box scheme). Let us consider vi > 0. After m wall 
reflections we find the particles in the box m. The true particle position is xi{£) = vit—Ana = X( j lee — Ana 
for even m = 2n (particle velocity ve) and xt(t) = 2a—{vzt—Ana) = 2a+Ana— xif ree for odd m = 2n+l 
(particle velocity — vg). Since the initial velocity distribution is symmetric and does not change in time, 
P{v, t) = Pq(v), we obtain for the distribution function for all v 

f(x,v,t) = e(l-^\ I jr f bee (x + 4na,v,t)+ f {lee (2a + Ana - x,v,t) ) . (10) 

^ / \n— — oo n— — oo / 
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The free Boltzmann distribution (O is symmetric under the transformation (—x, —v) —> (x, v). Therefore, 
we can simplify this expression to J4) 

f(x,v,t)=e(l-^) jr f boc (x + 2na,v,t) . (11) 

n— — oo 

We consider the Fourier series of the Boltzmann distribution with respect to the position coordinate and 
find for the Fourier coefficients (k = nir/a; n: integer) 

fk(v,t) = Yaf dx f(^ v ^ e " kX 

~ 2a 
1 

~~ 2a' 

The Boltzmann distribution function in the presence of hard walls and the free Boltzmann distribution 
function have the same Fourier coefficients at k = nir/a so that 



dxf i[eo {x,v,t)e- ikx (12) 

OO 

e~ lkvt P n (v) . 



/(*,«,*) = i-e(i-M) jr e i«w«)(— t) Po{v) _ (13) 



n— — oo 



With this result, we can calculate the probability distribution for the particle positions. We use the Fourier 
transformation for the initial velocity distribution, 

Po(v) = r %^P (V) , h(v) = [°° dve~ imi Po(v) (14) 

J-oo 27r J-oo 

with -Po(^) = Po(—v) an( l nn( l 

W(x,t) = ^<d(l-—] V e lxn7l/a [ dve- intV7,/a P (v) 







N 












a 






N 












a 



71— — OO 
OO 



e mxn ' a P {rmt/a) . (15) 

=— oo 

For each n^O, the argument of Po(t]) becomes large for t —> oo. We demand 

Po(t 7 ^oo) = 0. (16) 

This implies that the initial velocity distribution cannot contain a macroscopic number of particles with the 
same velocity. Then, for large times, the sum in ( fT5l ) reduces to its contribution from n = 0, 

w(x,t _> c») = ^-e f i - M") . (i7) 

2a \ ay 

For long times, the particles become equally distributed in the box with density po = N/V = N/ (2a), as 
expected from thermodynamic considerations. 

It is important to note that we know the particle positions and their velocities exactly at every instant 
in time. Despite this fact, the particle distribution becomes homogeneous for long times and 'relaxes' 
to its equilibrium distribution. For the particle positions this does not come as a surprise because there 
is no conservation law which prevents them to distribute equally in space. For the velocity distribution, 
however, momentum and energy conservation guarantee that it is independent in time for hard walls and our 
symmetric initial distribution. The relaxation of the velocity distribution is more involved, as we discuss 
in Sec.g] 
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3.3 Pressure on the hard walls 

The reflection of the particles at the right wall implies a momentum transfer AP to it. All the particles with 
velocity v > which reach the wall within the time span At transfer the momentum 2mv Q- Therefore, 

/•oo r-a 

AP = N dv(2mv) dxn(x,v,t) . (18) 

JO Ja-vAt 

For small At and in the thermodynamic limit, we find for the pressure in a one-dimensional box, p(t) = 
lim A t->o(AP)/(Ai), and with V = 2a 

f°° mv 2 — 
p(t)V = 2N dv — (2a)f(a,v,t) = 2N(E + H(t)), 

J -oo 1 

/oo 2 
Tnv 
dv—P (v) , (19) 



H(t) = / dv(mv 2 )e- m7TVt/a P (v) 



(20) 



The time-dependent terms, H (t), vanish for f ^ oo because the Fourier transformed velocity distribution 
Po(r]) and its derivatives vanish for large arguments. We find for t > 

H(t)=mf;(-1)-+ 1 ^'(^) , (21) 



n=l 



where Pq(i]) denotes the second derivative. 

The time-independent term provides the ideal-gas law (Q~|i in one dimension J2), p{t — > oo)V = 2NE 
where E is the average energy of a particle. For the Maxwell distribution, see eq. (0), we find E = k-gT/2 
so thatp(£ oo)V = Nk B T, eq. ©, results. 

3.4 Examples 

As examples, we consider an exponential velocity distribution and a step-like velocity distribution, 

P„(«) = ^e-l^ol , PM = T ^ R , (22) 

W«) = ^efi-M) , P ste p(.) = ^). (23) 

The particle density at the wall is obtained from (TT~5T > 

W(a,t) 



^ (^1 + 2 f^(-l) n P a (nnt / a^j . (24) 

With the help of Eq. (27.8.6) in Ref. and Mathematica (6) we find for t > (i tr = a/v Q ) 

°° (—1)" 1 
2aWU*,t) = 1 + 2^ ^L^ . ——— , ,25, 

»'.,m = i +2 f ( -ir25<^l=i-^. <26) 

n—l 1 ' 

Here, s(x) = x — 2n for x € [—1 + 2n, 1 + 2n] so that \s(x)\ < 1. The transient time i tr is given by 
the time a typical particle needs to travel through the system, t tr = a/vo- After this time, the density at 
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Fig. 1 Time evolution of the particle density at the wall, a, and of the pressure, b, for the step-like and exponential 
velocity distributions. Note that the average kinetic energy differs for the two distributions, 2E BtC p/ (rnvg) — 1/3 and 
2E exp /(mvi) = 2. 

the wall begins to approach its equilibrium value, NW(a, t —> oo) = po. For the exponential distribution, 
deviations are of order (t tl -/t) 2 for large times. For the step-like velocity distribution, we observe saw-tooth 
oscillations which die out for long times only proportional to (tu/i). This shows that the particles swap 
back and forth through the box and reach their homogeneous distribution only fairly slowly. 

The same analysis is readily carried out for the pressure. With the help of Eq. (27.8.6) in Ref. Q and 
Mathematica (6) we find for £ > {t tl - = a/vo) 



-Hstep (t) 



71=1 



1 - 3(nnt/t tl .) 2 
[1 + (n7rf/*tr) 2 ] 3 



UIVq 



3 + cosh(2£ tr /£) 



2[(£/£ tr )sinh(£ tr /£] 



mv. 



o 2 E(- 1 )" 



2cos(mr£/£tr) sin(n7r£/£ tr ) 



(nwt/tt 



mrt I £ t 



(n-Kt/ttrY' 



K£/£tr) , 3(s(£/£ tr )) 2 - 1 s(£/£ t r)(s(£/£tr) 2 - 1) 



3(£/£ t 



3(£/£ tr ) 3 



(27) 



(28) 



For large times, £ S> £tr, we find for the smooth exponential velocity distribution that H cx (t 3> £ tr ) ~ 
(£ tr /£) 4 , i.e., the pressure corrections vanish quickly as a function of time. For the step-distribution we 
observe again a slow decay with saw-tooth oscillations, H stcp (t 3> £ tr ) ~ (£/£tr) _1 - The density of 
particles at the wall and the pressure are shown in Fig.Q] 



4 Relaxation of particle velocities 

In order to relax the energy of the gas particles, we introduce scattering processes in which the gas particles 
exchange energy with the walls. 



4. 1 Dynamic wall and scattering events 

In order to simplify the analysis, we assume that the walls are perfectly symmetric. Each scattering event 
which occurs for the gas particle ge with velocity vt at the right wall (x = a) has a mirror event for the 
particle ~g e with velocity vj = —vt at the left wall (x = —a). In this way, our velocity distribution remains 
symmetric around v = at all times. In this section we are primarily interested in the relaxation of the 
velocity distribution of the gas particles. Therefore, we focus on the velocity, not on the position of the gas 
particles. 
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The right wall is made of L ^> N ^> 1 particles which move freely in the interval [a, a + A] (a 3> 
A > 0). The wall particles are scattered elastically at the interval boundaries. The boundary at x = a is 
transparent for the gas particles. The velocities of the wall particles are taken from a probability distribution 
D(u) = D(—u) which obeys the same restrictions as Pq(v), i.e., all of its moments exist and its Fourier 
transform, D{rj), decays to zero for large values, lim^oo D{rj) = 0. 

The gas particle gg (£ = 1, . . . N) which starts with velocity V£.o > at time t = undergoes a first 
scattering event at time t\ > with the wall particle with velocity Ui\. After the scattering process 
is completed, the gas particle has the velocity V£i < 0. At time t-i > t\ it is scattered off the left wall by 
the wall particle wi.2 whose velocity U£,2 is independent of the velocity of the wall particle W£ t i. The same 
applies to the subsequent scattering processes: The scattering partners of the gas particles are all different, 
we,n = wp.n' only if £ = I' and n — n'. In this way we assure that the wall has no memory: the outcome 
of a scattering event does not influence the distribution D(u). This assumption expresses the fact that the 
wall particles form a 'bath' for the gas particles. 

Let the wall particles have the mass M > m. Then, a scattering event between a gas particle with 
velocity v and a wall particle with velocity u results in 

v = /3u — av , 

u = rf3v + au , (29) 

m 1 — r 2 

r= M ' a = TT7 ' ^ = TT7 

with < r < 1, 0<a<l, and 1 < f} < 2. 



4.2 Master equation for the velocity distribution 

When we focus on scattering events at the right wall, we must distinguish three cases. We assume that the 
wall particle and the gas particle meet at the left wall boundary, x = a. 

- u < av / j3: 

The gas particle is reflected, v' — f3u — av < 0. 

- av//3 < u < v. 

After the first collision, the gas particle penetrates the wall because its velocity is still positive, u' > 
v 1 > 0. The wall particle is reflected at the wall boundary at x = a + A, and collides with the gas 
particle again. The gas particle then has the velocity 

v" = 7« - 2a(3u < , 7 = a 2 - r/3 2 , (30) 

and escapes from the wall region because v" < u" . 

- v < u: 

In this case, the wall particle escapes from the gas particle. The gas particle penetrates the wall and 
collides with the wall particle which was reflected at the wall boundary at a; = a + A. After the 
collision, the gas particle has the velocity v' = —/3u — av < 0. It cannot be reached by the wall 
particle again, v' < u'. 

If the average kinetic energy of the gas particles and the wall particles is comparable and in the limit 
m/M — > 0, the typical wall particle velocities are small compared to the typical gas particle velocities. 
Then, only the first case is relevant. 
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After all gas particles have been scattered n times, the probability distribution after the nth scattering 
process (n = 1, 2, . . .) is given by 



P J Jo P V P 



v'/(a+P) 



2a/3 V 2a P 



where we used the fact that our probability distributions are symmetric around the origin, P(n; v) = 
P(n\ —v), D(—u) = D{u). We can rewrite this equation in the form 

/>oo 

P{n;v)= dv'A(v, v')P(n - 1; v') (32) 



where the integral kernel A(v,v') is readily obtained from (l3TT l. Note that the wall has no memory. 
Therefore, the nth scattering distribution P(n; v) depends only on the distribution before the nth event, 
P{n — 1; v). This is the Markov property |3|. The probability distributions are normalized to unity. The 
condition /_ dvA(v',v) = 1 implies that no gas particle is lost which guarantees that the probability 
distribution remains normalized to unity, dv P(n; v) = 1. 

In order to elucidate the physical meaning of Eq. (l32l . we introduce the artificial time r = nAr 
(At = 1), and consider the limit of a large number of scattering events, n ^> 1, so that At t. 
We write 

P(n+l;v)-P(n;v)= dv' (A(v, v')P(n; v') - A(v', v)P(n; v)) . (33) 



In the limit n — > oo, the time r becomes a continuous mesoscopic time scale and we obtain a Master 
equation (3] for the velocity distribution P(t; v), 



8P(t;v) 



dv' (A(v,v')P(t;v')-A(v',v)P(t;v)) . (34) 



This is a gain-loss equation for the number of particles N(v) = NP(v) with velocity v. As has been shown 
in the literature |3j, the master equation (134-b has a unique stationary solution P sta t (v) = lim r ^oo P(t; v). 



0.4 



0.3 



> 
0^ 



0.2 



0.1 



— step-like velocity distribution 
--- exponential velocity distribution 



2 4 

v/u„ 



Fig. 2 Stationary velocity distribution for the gas par- 
ticles after an infinite number of collisions with the 
wall. The mass ratio between gas particles and wall par- 
ticles is m/M = 1/2. The velocity distribution of the 
wall particles is exponential and step-like, see Eq. i35i . 



In general, the stationary solution of the Master equation (f34b is not the Maxwell velocity distribution. 
As an example, we show the result for m/M = 1/2 for exponential and step-like velocity distributions of 
the wall particles in Fig. [2] 



\u/u \ 



2u 



D step (u) 



2uc] 



-9 1 



(35) 
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The stationary velocity distributions depend only on the choice of the velocity distributions for the wall 
particles, D(u), and not on the initial velocity distribution of the gas particles, Pq(v). The result is in- 
dependent of the initial conditions. Note that a Gaussian velocity distribution D(u) does not lead to a 
Maxwell velocity distribution P(v). 

The wall in the above example is not a 'physical' wall. According to the definition given in van Kam- 
pen's book |3|, a physical wall relaxes the gas particles' velocities to a Maxwell velocity distribution. 



4.3 Relaxation to a Maxwell distribution for a physical wall 

For a thermal reservoir, its properties must remain essentially unchanged by the interaction with the gas. 
If a wall particle transferred or received a substantial amount of energy in an individual scattering process 
it would carry information about the gas particles. Therefore, for a physical wall, we have to demand that 
the energy transfer in each collision is small compared to the typical kinetic energy of a wall particle. This 
can be guaranteed in the limit m/M — > 0. More specifically, we presume that the typical kinetic energy 
of a gas particle and of a wall particle are of the same order of magnitude mu t 2 yp ps MuL p . Then, the 
typical velocity of a wall particle is small compared to the typical velocity of a gas particle, i.e., u typ is of 
the order ^m/Mv typ <C Ut yp . 

In order to perform the limit m/M — > 0, we rewrite Eq. OTT l in the form 

, , , f°° D(x) fBx + v'\ f v ' /f} D(x) fv'-Bx 
P(v' > 0) = / dx— —P - + / dec— —P — 

J-v'0 a V a J Jv'/(a+(3) a \ a 

v ' a ^ D(x) f2aBx-v'\ 

dx-^-P — , (36) 



v'/(a+P) 



\l\ V 7 



where we used lim, woo P(n; v) = P(v). Now we are in the position to let m/M — > where a rj 
1 - 2m/M 4 1,^2(1- m/M) -> 2, 7 w (1 - 8m/ M) -> 1. D{x) is finite only in a region which 
is small compared to the typical values for v'. Therefore, we may safely extend the first integral over the 
whole real axis, and we may ignore the contributions from the second and third integral because D(v') = 
for typical values of v'. The excluded region, v 1 £ [0, (m / M)u t y P ], becomes vanishingly small in the limit 
m/M -> 0. 

The resulting integral equation can be simplified by Fourier transformation. We find, see eq. (fT4l . 

P{rf) = b{Bri)P{ar,) . (37) 

For the normalized, differentiable, and symmetric P{t)) and D(jj), the Taylor expansion of this equation 
gives 0(0) = 1, D'(0) =0,j8 = 2,a=l - 2m/ M) 

Pir,) = (p( V ) ^fP'iv)) + 2v 2 D"(0)P(v) ■ (38) 

Here, we made use of the fact that D(rf) is almost constant in the region where P(rj) is finite. This 
reflects our assumption that D(v) is a narrow distribution in comparison to P(v). The resulting differential 
equation is readily integrated. We use 

/>oo 

D"(0) = / duD{u){~m) 2 = -(u 2 ) , (39) 

J OO 

and define the wall temperature T via 

which corresponds to the equipartition theorem in one dimension. Note, however, that in our case k-oT/2 
simply is an abbreviation for the average kinetic energy of the wall particles. Then, we find from ( l38l 

P(r?) = e -'7 2fc BT/( 2m ) (41) 
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so that we finally obtain the Maxwell velocity distribution for the gas particles, 

™-™-yfjz5->(-*Sr)- <42) 

The average kinetic energies of the gas particles and the wall particles are the same, M(u 2 ) = m(v 2 ) = 
ksT, but their typical velocities differ by a factor yJrnjM, u typ = \J (u 2 ) = yJrajMvt yv with w ty p = 

yfW)> 

5 Numerical simulations for a deterministic system 

Finally, we present the results for a numerical simulation of gas particles in a box with a deterministic, 
physical wall. For simplicity we choose the left wall to be perfectly reflecting. 

5 . 1 Oscillator model 

We represent the right wall by harmonic oscillators of mass M with frequency ujq. Each gas particle gi 
(1= 1,2..., N) undergoes a sequence of collisions with the wall particles we.i (i — 1,2,..., Lf, Lt ^> 1) 
whose free motion is given by 

X£,i(t) = a - Acos(bj t + (f>£ ti ) , u*,i(i) = (Auj ) sin(w i + <Ae,i) ■ (43) 

The length A parameterizes the 'penetration depth' of the gas particle into the wall. We assume it to 
be a small fraction of the interatomic distance between the wall's atoms and, therefore, we work with 
A = 1 ■ 10~ 12 m = O.OlA. We choose a = 0.2m for the diameter of the box, a/ A = 2 ■ 10 11 . The typical 
phonon frequencies in a solid are ujq = 1.3- 10 13 Hz so that uq = Aojq = 13m/s. The phases S [0,2tt[ 
are chosen at random. This reflects the fact that the motion of the wall particles is uncorrected. We shall 
comment on the case of equal phases in Sect. 15.21 

For long enough times, we expect that we recover the situation of Sect. 14.31 The velocity distribution of 
the wall particles is given by (To = 2it/ujq) 

1 f T ° 11 
= — dtS (u - u sm{uj t + 0)) = for \u\ < u . (44) 

The Fourier transformed velocity distribution is given by a Bessel function D(rj) = Jo{uorj) which decays 
to zero for large arguments, D(rj 1/uo) oc 1/y/fjuo. 

From d40b it follows that the wall temperature is given by /cbT = Muq/2 because (u 2 ) = Uq/2 for 
the velocity distribution (l44l i. Therefore, at room temperature we find M « 29 • 10 3 u (u: atomic mass 
unit), i.e., of the order of O(10 2 ) atoms in the wall constitute a 'wall particle'. As our mass ratio we set 
m/M = 10~ 3 , i.e., the gas particles have the mass m = 29 u (nitrogen molecule, N2). 

We expect that the stationary velocity distribution of our gas particles becomes a Gaussian distribu- 
tion d42b with variance (v 2 ) = k B T/m = Mj{2ra)u\ = 500uq. We thus expect v typ = yj (v 2 ) = 23u ps 
3 • 10 2 m/s which agrees with the typical velocity of gas particles in a container at room temperature. The 
time for a typical particle to perform a closed loop through the box is it yp = 4a/i>t yp ps 1.3 • 10~ 3 s. A 
good estimate of the number of scatterings C(t) which a gas particle has undergone after the time t is given 
byC(t) = O(t/(10ttr)). 

At time t = 0, we distribute our A^/2 = 1 • 10 7 gas particles equidistantly in the velocity interval ]0, vq], 
Pq(v > 0) = 9(1 — v/vo)/(2vo), with vq = lOito = 130m/s. The results for negative velocities are 
obtained by a reflection at the axis v = 0. Initially, all gas particles start at x = 0, f(x,v,to = 0) = 
P (v)5(x). 
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Fig. 3 Velocity distribution P(v > 0, t) in the oscillator model for N/2 = 1 ■ 10 7 particles for v > at times to = 0, 
ti = 0.1s, t 2 = Is, t s = 10s, corresponding to a number of C(t ) = 0, C(ti) = ©(lO 1 ), C(t 2 ) = O(10 2 ), 
C(t 3 ) = O(10 3 ) wall scatterings. Full lines: Maxwell distribution with (v 2 ) = Mul/{2m). 



5.2 Results 

5.2.1 Random oscillator phase 

In Fig.[3]we show the velocity distributions P(v > 0,t) at times to = 0, t\ = 0.1s, = Is, = 10s, 
corresponding to a number of C(t ) = 0, C{t{) = ©(lO 1 ), C(t 2 ) = O(10 2 ), C(t 3 ) = O(10 3 ) wall 
scatterings. We see that the velocity distribution continuously transforms from the initial step-distribution 
to the desired Maxwell distribution. The calculated variance of the essentially stationary distribution at 
t = ts is a 2 = 502mq, the expected value is (v 2 ) = Muq/ (2m) = 500uq = k^T/m. Apparently, our 
oscillator model fulfills the criterion for a physical wall because it leads to a relaxation of the initial velocity 
distribution to a Gaussian. For very small velocities, we observe minor deviations from the Gaussian shape. 
This can be attributed to the fact that our ratio m/M is finite so that the condition |u| > (m/M)uo is 
violated for very small \v\, cf. Sect. 14.31 

We started from other non-singular velocity distributions with -Po(M > vo) = 0, and also recovered 
the Maxwell distribution as stationary state of our simulation. In addition, we verified numerically that the 
Boltzmann distribution for long times t is given by 

f(x,v,t» i tr ) = ^ ( x - ^) /m(") • (45) 

As expected from our discussion in Sect. [3] the gas particles are equally distributed in the box, and there is 
no correlation between their position and their velocity. 
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5.2.2 Fixed oscillator phases 

Finally, we address the issue of whether or not our model can be constrained further. In our simulation 
results shown in Fig. [3] the phase (j>i^ of each wall oscillator is chosen at random. If this phase is fixed, 
e.g., to (j>t i = 0, all oscillators start with the same phase at t = 0. When the gas particle gi reaches the 
wall, vi$tt\ — a — A, the first oscillator has the contact phase 7^1 = (t^.i mod To) (To = 2tt/luq). After 
its scattering(s) off the wall particle, the gas particle leaves the wall region at xt t i with velocity —V£ t i < 0, 
is reflected at the left wall, and returns to the right wall boundary at time tip = (2^1 + 3a — A)/v£.i 
with velocity vi : \. The contact phase of the second oscillator is 70 = (te,2 mod To). In general, if tt n 
is the time at which the gas particle reaches the boundary of the right wall for the nth time tg !n = 
(x£. n -i + 3a — A)/vi t n-i with velocity v^„_i > 0, the contact phase is given by 7^„ = {t^ n mod To). 
The scattering off the wall can be seen as a dynamical map (n > 1), 

$ : («n,7») !->• («n+l,7n+l) ( 46 ) 

for the velocity of the gas particles and their contact phase. The map $ belongs to the class of Fermi-Ulam 
accelerators (7). 
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Fig. 4 Stationary velocity distribution for the gas particles after a large number of collisions with the dynamic wall. 
The maximal wall particle velocity is tio = 13 m/s. The mass ratio between gas and wall particles is m/M — 10 -3 , 
the size of the box is a = 100A, and the particles evolve according to the deterministic map <E>. The largest peak is a 
fix point of the map, v p = i>+(3) = 277 m/s. 

In general, the Fermi-Ulam accelerator map does not lead to a Maxwell distribution for the gas particle 
velocities. In order to see this, we note that $ has stationary points, 7n+i) = (v n> 7n)- They follow 

from the solution of the equations 

v n = av n - /3u(t n ) , 

2(a + x(t n )) 2nR 
t n+ i-t n = 1 ^ Kn,} = (R>1: integer) , (47) 

where u(t n ) = Aluq sin(woin) and x(t n ) = a — A cos(ujot n ). One finds for R > R m j n = 2ra/ (Aw) ^> r 
that v±(R) = uo[2a/A± y/l — (2ra/ (AttR)) 2 ]/(ttR). Moreover, the map supports stationary loops with 
two and more elements. As we have tested numerically for a = 100A, some of the fix-points and loops are 
attractive with a large basin of attraction around them (e.g., R = 3). Thus, the stationary distribution P(v) 
depends on the initial distribution Pq(v) and consists of a few peaks only, see Fig. [4] In this simulation, we 
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also observed that these peaks belong to small values of R. The stationary solutions for large R, R 1, 
do not generate attractive fix-points. 

The situation changes drastically for large raj A, i.e. when we choose i? m ; n = 0(ra/A) = 10 8 . As 
expected for a Fermi-Ulam accelerator model and confirmed by our numerical investigations, the map $ 
in d46b is chaotic for ra ^> A; see Refs. (8) for a recent introduction to dynamical systems and chaos. 
For example, tiny numerical errors make it impossible to retrace the gas particles' trajectories when we 
reverse the time evolution after a large number of wall scatterings. This implies that the deterministic map 
in Eq. d46b is acceptable for the thermalization of the gas particles in our toy model as long as we demand 
that the box is (mesoscopically) large compared to the penetration depth of the gas particle into the wall, 
ra 3> A. A more detailed investigation of the map $ is beyond the intentions of this work. 

6 Summary and conclusions 

In this work, we studied non-interacting point particles in a box. Since we assume that the initial conditions 
factorize in the three spatial dimensions, we restricted ourselves to a one-dimensional setup. As initial 
condition, all gas particles start in the center of the box with a non-singular velocity distribution. The 
particle positions 'thermalize' after a large number of scatterings at the box walls, i.e., the gas particles 
become equally distributed in the box. This is due to the fact that there is no conservation law for the 
position coordinate. The notion of equal spatial distribution is valid on mesoscopic length and velocity 
scales Ax — a/ y/N, Au = i>t yp /\/iV, which requires a large number of gas particles, N ^> 1, 

The relaxation of the gas particles' velocity distribution to a Maxwell distribution is more subtle because 
of energy and momentum conservation. We model our walls as a large collection of independent harmonic 
oscillators, L 3> N . Each gas particle encounters a harmonic oscillator only once (absence of a wall 
memory). The Markov property is an essential requirement for finding a stationary velocity distribution 
for the gas particles. In order to find a Maxwell velocity distribution, we must demand in addition that 
there is a negligible energy transfer between the gas and the wall particles in each collision; otherwise, 
a single event could leave a trace of information of the gas particles in the reservoir. The condition of 
small energy transfers leads to the requirement that the wall particle mass M is large compared to the gas 
particle mass m, and that their typical velocity ut yp is small compared with a typical gas particle velocity 
Utyp. Their average kinetic energies must be the same, M (u 2 ) — m(v 2 ). Finally, the oscillation amplitude 
of the wall particles A must be small compared to the size of the box a so that the map of the Fermi- 
Ulam accelerator guarantees a chaotic dynamics. When these 'thermodynamic' conditions are fulfilled, the 
velocity distribution of the gas particles approaches a Maxwell distribution for large times. The evolution 
of the full system of gas and wall particles is fully deterministic but the result, a Boltzmann distribution for 
the gas-particle subsystem, is naturally obtained from statistical considerations. The 'temperature' of the 
reservoir is defined by the average kinetic energy of the wall particles. M(u 2 ) = k&T. 

Apparently, our toy model does not describe the details of the relaxation process properly. For example, 
the scattering of gas particles at the walls must be described quantum mechanically. As a consequence, 
not every scattering process is inelastic. In fact, most of them are elastic whereby the probability for 
an elastic reflection is determined by the Debye- Waller factor. Therefore, the relaxation times could be 
(much) longer for non-interacting gas particles than given here. On the other hand, the interaction between 
gas particles and their motion through a three-dimensional vessel reduces the relaxation time drastically 
because the inter-particle scattering opens paths for the exchange of energy and momentum which are not 
included in our toy model. In any case, the relaxation time also depends on the shape of the initial velocity 
distribution. If it is peaked around some velocities, it takes much longer to equilibrate than for the step-like 
distribution used in our example. 

Despite its numerous shortcomings we hope that our work will serve its intended main purpose, namely, 
to provide a simple case study for beginners in statistical mechanics. Our toy model illustrates that the 
results of measurements on subsystems of thermodynamically large deterministic systems can equally be 
obtained from thermodynamic considerations. 
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